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<^ ■ Abstract 

^ . The LS'-sequences of points recently introduced by the author are a gen- 

I eralization of van der Corput sequences. They were constructed by reorder- 

CL|' ing the points of the corresponding L5-sequences of partitions. Here we 

■^ . present another algorithm which coincides with the classical one for van der 

QQ [ Corput sequences and is simpler to compute than the original construction. 

This algorithm is based on the representation of natural numbers in base 

L + S and gives the van der Corput sequence in base b it L = b and 5 = 0. 

f-H , In this construction, as well as in the van der Corput one, it is essential the 

^ I inversion of digits of the representation in base L + S: in this paper we also 

fH ■ give a nice geometrical explanation of this "magical" operation. 
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g : 1 Introduction 



This paper contributes to enhance the connection between the classical theory of 



O ■ uniformly distributed sequences of points in [0, 1 [ and the related area of uniformly 



distributed sequences of partitions of [0, 1 [, by showing how to reorder the points 

of the LS-sequences of partitions in order to get the LS-sequences of points, both 

introduced by the author in [5]. If L = b and 5 = 0, the former generalize the 

rN ■ sequences of partitions in base b, the latter generalize van der Corput sequences 



c^ ■ in base b. 



The construction presented here is based on the representation of positive in- 
teger numbers in base L-\-S and coincides with the construction of van der Corput 
sequence in base b in which, as it is well known, the representation of positive 
integer numbers in base b is essential, as well as the inversion of these digits and 
the corresponding radical-inverse function 0^. The most interesting contribution 
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of this paper consists in the geometrical explanation of the above mentioned in- 
version of digits, which appears not only for the van der Corput sequences, but 
(together with some more specific features) also for all the LS-sequences of points. 
Let us begin by giving some definitions (for a complete overview on uniform 
distribution and discrepancy see [20] ). 

Definition 1.1. Given a sequence {Qn} of finite subsets of [0, 1[, with 2„ = 
{jj , . . . , j}" }, we say that it is uniformly distributed (u.d.) if 

limi£/(3;W)= f'f{x)dx (1) 

for every continuous (or Riemann integrable) function / defined on [0, 1 [. 

Definition 1.2. Given a sequence {jc„} of points in [0, 1 [, we say that it is uniformly 
distributed (u.d.) if the sequence of sets 2„ = {^i, . . . ,x„} is uniformly distributed. 

Definition 1.3. Given a sequence {;r„} of finite partitions of [0, 1[, with ;r„ = 
{ bz-i 5 yi [,!<''< tfi}, Jo = and y}" = 1 , we say that it is uniformly distributed 
(u.d.) if the sequence of sets Q^ = {y", . . • , j}„ } is uniformly distributed. 

These definitions can be easily extended to higher dimensions, but the results 
of this paper concern only the one-dimensional case. 

Sequences of partitions and sets have been considered rather early in the theory 
of uniform distribution. Partitions in base 2 have been considered in 1935 by van 
der Corput in [25 1. Knapowski [19| studied in 1957 the sequence of partitions of 
[0, 1 [ into n equal parts and its generalizations. Hammersley introduced van der 
Corput sequences in any base b and introduced a particular sequence of subsets of 
R" which proved to be efficient when evaluating multidimensional integrals with 
Monte Carlo type methods (see [|T2l for these applications). The L5-sequences in 
the multidimensional case have been studied in [|61 and [|3|. 

More precisely, the sequence of partitions in base b, with b E'N, b > 2, is 

defined by 

"/— 1 i 



b" ' b" 



l<i<b"\ (2) 



while the corresponding sequence of points, called the van der Corput sequence in 
base b, is obtained reordering in an appropriate way the points determining them. 
We will not discuss this aspect in detail now, as this procedure turns out to be a 



special case of the algorithm we will propose in last section of this paper. We only 
mention Faure's scrambled van der Corput sequences ( lfT3][T4]| ') and the previous 
extensions to higher dimension by Hammersley [|T6l and Halton ifTSl . 
Knapowski considered sequences of points 

12 ai - 1 1 2 a2 - 1 1 a„ - 1 "1 

, , • ■ ■ 7 7 7 • • • 7 7 7 • • ' 7 7 ' ' ' f \^ } 

a\ a\ a\ a2 a2 aj an an ) 

related to equipartitions of [0, 1 [. He analyzed the conditions under which the 
sequence ([3]) is u.d. (see also p23]). 

Hammersley considered sets of N points in R"^ of the kind {x" , ^Cj 5 • • • ^d ) 

(n) (n) 

for n <N, where x\ = fj, while for 2 < ^ < J, x^ is the n-th term of the van der 
Corput sequence of points in base b^^i, the {k— l)-th prime number. 

From the definitions of uniform distribution of sequences of sets, points and 
partitions it is clear that, when we want to evaluate the integral of /, the sequences 
of points provide a more flexible tool, as we can choose in advance any number A'^ 
of points Xn and afterwards, if we want to improve the approximation, we can add 
any number of additional points using all the previously calculated values of /. 

On the other hand, if we want to increase the number of points and we are 
dealing with a sequence of sets (or partitions), there are two possibilities. The 
worst case is when Qn and Qn+k have few points in common. In this case all the 
previously calculated values of / have to be thrown away. The situation is better 
if Qn C 2/7+1 . But even this case has a drawback, as we are limited in the choice 
of the number of points in which to evaluate the function /, as the averages ([I]) 
have no meaning for values between r„ and tn+i- 

From this observation stems the main motivation of this paper: given a se- 
quence of u.d. sequences of partitions, it is highly desirable to associate to it a u.d. 
sequence of points having the lowest discrepancy. In [5| the author introduced a 
class of sequences of partitions (the so-called L5-sequences), which includes the 
sequence of partitions in base b defined in ([21). In the same article she proposed an 
algorithm which gives a uniformly distributed sequence of points reordering the 
points determining the corresponding sequence of partitions. When the sequence 
of partitions has low discrepancy, the discrepancy of the corresponding sequence 
of points is the same up to a constant times a logarithmic therm (as we will recall 
in the next section), and this is optimal too. If L > 2 and 5 = 0, the sequences 
reduce to the van der Corput sequences in base L. 

The purpose of this paper is to present another algorithm which constructs the 
same sequences but in a more efficient way. 



2 Preliminaries 

Kakutani introduced in the seventies the sequence of successive a-refinements of 
a partition of the unit interval [[T8l . 

Definition 2.1. Given a finite partition K of [0, 1[ and given a G ]0, 1[, its a- 
refinement, denoted by an, is the partition obtained by subdividing all the inter- 
vals of K having maximal length in proportion a and 1 — a. We denote by a";r 
the ot-refinement of a"^^;r, with a^Tt = an. 

He proved the following result. 

Theorem 2.2. If 7t = (O = {[0,l[} is the trivial partition of [0, 1 [, the sequence 
{a"n} is uniformly distributed. 

Kakutani's result got a considerable attention in the late seventies and early 
eighties, when other authors provided different proofs of Kakutani's theorem [HI 
and also proved several stochastic versions, in which the intervals of maximal 
length are split according to certain probability distributions ( fl2T]| . [|22l . ||26l , and 

a). 

The paper [jTOl extended the notion of uniform distribution of a sequence of 
partitions to probability measures on complete separable metric spaces. 

In [8] Kakutani's splitting procedure has been extended to the n-dimensional 
cube with a construction which is intrinsically higher-dimensional, and in IfTTI it 
has been extended to fractals. 

In im the authors present a von Neumann-type theorem for sequences of par- 
titions, showing that a sequence of partitions having an infinitesimal diameter can 
be reordered to a sequence of partitions which is uniformly distributed. 

The paper [|27l introduced the following generalization of the splitting proce- 
dure. 

Definition 2.3. For any non trivial finite partition p of [0, 1 [, the p -refinement of 
a partition TZ of [0, 1 [ (denoted by pTZ) is obtained by subdividing all the intervals 
of K having maximal length positively (or directly) homo thetic ally to p. If for 
any n G N we denote by p";r the p-refinement of p"^^;r, we get a sequence of 
partitions {p";r}, called the sequence of successive p-refinements of n. 

Obviously, if p = {[0, a], [a, 1 [}, then the p-refinement is just Kakutani's a- 
refinement. 

As in Kakutani's case, we can iterate the splitting procedure. The p-refinement 
of p;r will be denoted by p^;r, and the meaning of p";r, for n G W, is clear. The 
following result generalizes Theorem 12.21 



Theorem 2.4. The sequence {p"(U} is uniformly distributed. 

The recent paper by Aistleitner and Hofer [|2] provides a complete solution 
to the question, posed in [|27ll . for which partitions K the sequence of partitions 
{p";r} inu.d.. 

The following result in [|27l connects u.d. sequences of partitions to u.d. se- 
quences of points. 

Theorem 2.5. Given a u.d. sequence of partitions {7r„}, a random, sequential 
reordering of the points determining the partitions is a u.d. sequence of points 
with probability one. 

A sequential ordering of the sequence of sets Qn = {jj , . . . ,y}' } is the se- 
quence whose first ti points belong to Qi, the next ?2 points belong to Q2, and so 
on. 

The limit of this theorem is that it does not provide an explicit algorithm to 
construct a u.d. sequence of points associated to {;r„}. As we will see in the next 
section, for the family of L5-sequences introduced in [5], which are a subfamily 
of r/?o-refinements, this algorithm is explicitly given. 

We finish this section by giving some very important definitions. 

Definition 2.6. Given a sequence {Qn} of finite subsets of [0, 1[, with Qn — 
{y" ,. . . , yt }, its discrepancy is defined by the sequence 



D{Qn) = sup 

0<a<b< 1 



tn j^i 



If we look at Definition 1.1, we see that discrepancy gives us the speed with 
which the averages f-E/Li X[a,b[{yi ) converge to the integral of / = X[a,b[- Here, 
as usual, X[a.b[i') denotes the characteristic function of the interval [a,b[. 

Definition 2.7. Given a sequence {Qn} of finite subsets of [0, 1[, with Qn = 

(n) (n) 

{y\ , . . . , yt„ }, we say that it is has low discrepancy if there exists a constant 
C such that tn D{Qn) < C. 

It is well known that the upper bound cannot be improved and that it is optimal 
(see [20], for instance): the Knapowski sequence of sets has low discrepancy. 

However, if the sequence of sets Qn is generated by a sequence of points as in 
Definition 1.2, it is well-known that the corresponding discrepancy cannot be so 
fast (see [12411 ). 



Definition 2.8. Given a sequence of points {xn} of [0, 1 [ and 2„ = {^i , . . . ,a:„}, we 
say that it is has low discrepancy if there exists a constant C such that nD{Q„) < 
Clogn. 

This bound is optimal and it is achieved, for example, by the van der Corput 
sequence in base b. 

Let us mention that a further contribution to the theory of p -refinement has 
been given by Drmota and Infusino [11] who provided estimates of the discrep- 
ancy of sequences of p -refinements and, in particular, of L5-sequences of parti- 
tions. 



3 L5-sequences 

In this section we recall the definition of LS-sequence of partitions {p'lg] and 
LS-sequence of points {(^^5} introduced in [|5l . 

Definition 3.1. Let us fix two positive integers L > 1 and S > with L-l-5 > 2 
and let 7 be the positive solution of the quadratic equation Lx -\- Sx^ = 1 (if S = 
the equation is linear). Denote by pis the partition defined by L "long" intervals 
having length 7 followed by S "short" intervals having length 7^ (if 5 = all the 
L intervals have the same length 7 = 1/L). The sequence of the successive Pl,5- 
refinements of the trivial partition O) is called LS-sequence of partitions and is 
denoted by {p£ ^ co} (or {p£' 5} for short). 

Note that for L = b and 5 = we get the sequence of partitions in base b 
defined by . 

If we denote with ?„ the total number of intervals of p£ ^, with l„ the number of 
long intervals and with Sn the number of short intervals, it is very simple to see that 

tn — In+Sn, In = Lln-l + Sn-l and Sn = 5/„_l. We dcduCC that tn = Ltn-\ +Stn-2 

(with to = I and ti = L + S) which has the explicit solution 



1+57 /ly Sy-Sr^ 



-57)". 



" 1+572 Vr/ l+Sy^ 

We note that /„ and 5„ satisfy the same difference equation, with initial conditions 
lo= l,li = L and sq = 0, si = S, respectively. 

The following estimates from above and from below of the discrepancy of 
L5-sequences of partitions have been given in [|5]- 



Theorem 3.2. i)IfS<L there exist ci,c\ > such that for any n G N 

c\<tnD{pls)<Cl. 

ii) IfS = L+l there exist C2, C2 > such that for any n G N 

C2i0gtn < tnD{pls) < C2logf„. 

Hi) IfS>L + 2 there exist C3, C3 > such that for any n G N 

cW <tnD{pls)<C3t'„-\ 

Where l-T = -'^>0. 

We emphasize that each L5-sequence of partitions with S <L has low discrep- 
ancy. 

Let us now recall the definition given in [5] of the sequence of points {^1^} 
associated to the sequence of partitions {p'lg}- 

Definition 3.3. We fix a sequence of partitions {plg}- For each n > 1 we define 
the two families of functions 

(pj"+^\x)=x + if+^ and (pl^j'^\x)=x + Lf+^+jf+^ (4) 

for 1 < / < L and 1 < j <S —I (if 5 = we have only the functions (p^ ' ). We 
denote by A| ^ the set made by the t\ left endpoints of the intervals of p^ ^ ordered 

by magnitude. If A2 ^ = ( ^l" ,. .., ^i^ ) denotes the set of the r„ points defining 

p£ g, ordered in the appropriate way, the set A2^^ , which consists of a reordering 

of the points defining p'^'^^ , is obtained as follows: 






,pr"(i;"').....'pr"«,i"'),...,?r"(i;"'),...,?r"(«,::'). 
<,+''(e')....,^&'^ir)'---.'pSA'(ir).---.<';--'i(i,!:'))- <« 

The sequence obtained by reordering successively the points of pjj^^ \Pls ^^ 
called LS-sequence of points and is denoted by {^l^}- 



Let us observe that if L = b and 5 = 0, the sequence reduces to the van der 
Corput sequence in base b. 

In [5J we provided the following bounds for the discrepancy of these se- 
quences of points. 

Theorem 3.4. i)IfS<L there exists ki>0 such that for any N eN 

ND[^ls^^ls^---^^L,s)<kilogN. 
ii) IfS = L+\ there exists ^2, Cj > such that for any N eN 

c'^XogN < NdI^^Is. ^Is^ ■ ■ • , &) < hWN. 
Hi) IfS>L + 2 there exists ki,,C2 > such that for any N eN 

c',N'-'<Nd[^Is.^Is.---.^l,s) < k3N'-'logN, 

wherel-T = -'^>0. 

Theorem 13 .2 1 and Theorem 13.41 show that if the L5- sequence of partitions has 
low discrepancy, the corresponding L5-sequence of points has low discrepancy, 
too. Moreover, Cj and c'^ denote the same constants in both theorem as, of course, 

whenever A^ = f,j for some n, it follows thatDK^^,!^^^, ... ,(^^^j = D{plg). 

Note that the lower bounds follow from Theorem 13.21 In fact we strongly 
believe that the discrepancy of all the L5-sequences of points coincides with the 
upper bound given by Theorem 13.41 

4 Technical results 

In this section we present an algorithm to construct L5-sequences of points {^ig}, 
which is formally independent from the concept of sequence of partitions: we do 
not need to reorder the points of each partition, as it is shown in Theorem 15. 31 and 
therefore it is simpler to use then the algorithm introduced in [5]. Of course, the 
connection between each {^'Ig} and the corresponding {p'lg} lies behind all the 
procedure, and this geometric interpretation represents one of the most interesting 
aspects of this new family of sequences. 



In the case L = b and 5 = 0, this algorithm will produce the van der Corput 
sequence in base b. 

From now on we fix L, 5 G N and 7 such that Ly+Sy^ = 1 . 

In order to introduce this algorithm, we represent the points of the sequence 
{^Ls} writing the elements of each set A'l g defined by (O as images of the point 
X = under the compositions of suitable functions depending on L and S we are 
going to define. 

For this purpose we introduce the functions i//, as follows: for every < ? < 
L — 1 we consider 

V^; W = yx + iy restricted to < x < 1, (6) 

while for every L< i < L + S— I and 5 > 1 we define 

Vi (x) = yx + Ly+{i-L)y^ restricted to <x<y. (7) 

If 5 = 0, we will consider only the functions ©. 

We observe that the functions (|6l) map [0, 1[ onto [iy, (/+ 1)7[ and the func- 
tions dV]) map [0,7[onto [Ly+{i—L)'f-,Ly+{i—L+\)y^[. Therefore, the compo- 
sitions I/// o x^fj are not defined if and only if L<i< L+S— I and 1 <j<L-|-5— 1. 

Let us denote by E^s the set consisting of all the pairs of indices which corre- 
spond to the "forbidden" compositions, i.e. 

EL,s = {m+l,...,L + S-l}x{l,...,L + S-\}. (8) 

If 5 = 0, the first factor is empty, so E^s = ©• 

In order to present our main result, we have to give the explicit expression of 
the compositions V^/i/j...;,, = ¥h °Wi2° ■■■ Win of the functions Q and ([7]). 

Lemma 4.1. For any n G N and any n-tuple {i\i2...in) of elements of the set 
{0, 1, . . .,L + 5— 1} such that (ihJh+i) ^ E^s for any \ <h <n — l, we have 

¥hi2...M=r^+t^k"/': (9) 

fe=l 

where b^ = ik if ^h ^ L—1, while bk = L+{ik~L)y if L <ik <L + S —1. 

Proof. We prove ^ by induction on n. Of course, we need to consider only the 
case 5 > 1. 

If n = 2, because of ^ and (|7]) we have only three kinds of compositions 

¥h,l2- 



a) if < i'l , /2 < ^ - 1, we have yA,; ^i^ (x) = 'fx + LLi ^^ f^ where hk = 4; 

b) if < ?'i < L — 1 and L</2<L + 5 — l,a simple calculation gives dH), with 
h\ =i\ and b2 = L+{i2 — L) 7; 

c) if L < /,■ < L + 5 — 1 and < i2< L—l,we obtain the same expression, where 
bi = L+ (z'l — L)7and b2 = 12- 

Suppose now that (fTTI) holds for n > 1 and let us prove it for n+ I, namely 
that 

¥n^2-,,inA^) = f^'^ + "t^ky' (10) 

k=l 

for all in+i 7^ 0. Of course, at this stage we have to distinguish between the two 

kinds of functions Q and (|7]). 

If < in+i < L — 1, we simply obtain (flOl) with bn+i = in+i- 

If L < i^_|_i < L + 5— 1, simple calculations give again (fTOl) . with Zj^+i = L + 

('«+i ~^)7- So the induction is complete. D 

Remark 4.2. We point out that ifS = 0, the coefficients b^ in formula dP]) are all 
equal to i^ as Eis = 0- 

As a consequence of the previous lemma we obtain the following result. 

Proposition 4.3. The first t^ points of the LS -sequence {(^£5} are 

AZ,5 = {V^n.2...J0)}, (11) 

where {ihJh+i) ^ E^s for any 1 <h<n — \ and {inJn-iT ■ -ih) are the n-tuples 
of elements oftheset{Q,\,...,L + S—l}, ordered with respect to the magnitude 
of the numbers inin- \- ■ -h in base L + S. 

Proof. Taking ^ and (|7]) into account, we see that 

Ai,5 = {V^o(0),V/i(0),...,v/L+5-i(0)}, 

which are exactly the L + S left endpoints of the intervals of pi g ordered by mag- 
nitude, so (fTTI) is true for n = 1. 

Fix any integer n> I and suppose (fTTI) is true for n. 

If in+i = 0, as V/b(0) = we have 

and we simply observe that this way we get all the points of A^g. All the new 



points of A"Y \ A? c are obtained when i^+i 7^ as follows. 
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If 1 < ?„+i < L we have 

n 

according to dU). 

If L+ 1 < in+i < L + S— 1, we set jn+i = in+i —L, with 1 < j„+i < 5— 1, and 
write 

k=l 

= Wni2...d0) +Lf+'+jn+if+' = 'pi;;;;(V^u.,... JO)). 

At this point we recall that the set A^"^^ \ A^ ^ is made by (L + S — 1 ) /„ points 
as we evaluate the L + S — I functions defined by (Hj) at the first /„ points of A'l ^ 
(see ©). 

It remains only to prove that the total number of the points of A^^^ \ A2 5 we 

add in the way described above is exactly {L + S — 1)/,,. We denote by d^^ the 
total number of n-tuples containing two consecutive digits of the set Eis defined 

by ©. Of course, we have 4^ = (L + S)" - 1^. 

We observe, in fact, that the points of A^^^ \ -'^l 5 come from L + 5 — 1 blocks 
of all the {n + 1) -tuples of elements of the set {0, 1, ... ,L + 5 — 1}, which are 
(L + S)", ordered by magnitude of the numbers in base L + S. Moreover, the 
total number of not admissible {n + 1) -tuples which are contained in these blocks 
is given by the difference between {L + S — l)(L + 5)" and the total number of 
forbidden ones, which is exactly ai^ — j]-"^. In other words, 

(L + 5-l)(L + 5)"-(ji;;^)-jW] 

= {L+s-\){L+sr-{{L+sr^'-t„+i-{L+s-ir+tn) 

— tn+l -~hi = {L + S — \)ln- 

The proposition is completely proved. D 
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5 The new construction 

We are now ready to give the construction of L5-sequences of points, which is 
deeply connected to the explicit geometric representation given in Proposition 
|4.3[ More precisely, we want to show how to get the points of these sequences 
starting from the representation of natural numbers in base L + S. This algorithm 
is very simple and allows us to compute directly the points of each {^[g}- 

This algorithm is, in fact, a generalization of the algorithm which produces the 
van der Corput sequences. 

Any natural number n> I can be expressed in base i» > 1 as 

M 

n=Y,ak{n)b\ (12) 

k=Q 

with ak{n) E {0, 1 , . . . , Z? - 1 } for all < A: < M, with M = [log^ n\ (here and in 
the sequel [ ■ J denotes the integer part). 

The expression (fT2l) leads to the representation in base b 

[n]h = aM{n)aM^i{n)...aoin). (13) 

If n = 0, we write [0]^ = 0. 

The representation of n in base b given by (fT2)) is used to define the radical- 
inverse function (j)h on N which associates to the string of digits (fT3] ) the number 

M 

(l>b{n) = Y,akin)b-'-\ (14) 

whose binary representation is 0.ao{n)ai{n) . . .aM{n). 

Of course < ^i,(«) < 1 for all n > and the b-radix notation of 0i,(n) actually 
is [^b{'^)]b = 0.ao{n)ai{n) .. .aM{n)- 

The sequence {<^fo(«)},j>o is the van der Corput sequence in base b. 

Definition 5.1. We denote by N/^ 5 the set of all positive integers n, ordered by 
magnitude, with [n]/^+5 = aM{n)aM-\{n) ...ao{n) such that {ak{n),ak+i{n)) ^ 
El^s forallO<yt<M-l. If 5 = 0, we have Nl,5 = N. 

Definition 5.2. For all n E Nl^s we define the LS-radical inverse function as fol- 
lows: 

M 
k=0 

12 



where dk{n) = ak{n) if < ak{n) <L—l and dk{n) = L+ 7(a^(n) — L) if L < 
ak{n) <L + S—l. (If 5 = 0, (fT5l) coincides with the radical inverse function (IT41).') 

Taking the above definitions into account, we are able to present the main 
result of this paper. 

Theorem 5.3. Any LS-sequence of points coincides with the sequence {(^l,s{^)} 
defined on 'H^s- 

Proof. Fix n G Nl,s- Taking Lemma 1431 and Definition 15 . 21 into account, we write 

which proves the theorem. D 

Remark 5.4. If L = b and 5 = 0, we have ^t.o = ^b <^nd, therefore, the b, 0- 
sequence is the van der Corput sequence in base b. All the results contained in 
this section can be re-written in this particular case, a. 

Example 5.5. In this example we discuss the 1, 1 -sequence of points in order to 
better illustrate the algorithm described in this section and to underline analogies 
and differences with the van der Corput sequence in base 2. For yet another appo- 
rach to the 1, 1-sequence of points, see fJ\. 

We consider the Kakutani-Fibonacci sequence of points {(^f i) corresponding 
toL = 5=l and7=^(V5-l). See Fig. 1. 

1 6 4 3 8 2 7 5 

I 1 1 1 1 1 1 1 1 

y^ y3 y2 y^ -U y^ T y _j_ y+ y _j_ y3 1 

Figure 1: The first 8 points of {(^"j}. 

The functions defined by ([8]) and ^ reduce to \]/o{x) = yx for < jc < 1 and 
Vi (x) = 7.x: + 7 for < .x: < 7 See the Fig. 2. 

According to Definition 15. 2[ Ni 1 is the set of all natural numbers n such that 
the binary representation ([131 ) does not contain two consecutive digits equal to 1. 
Moreover, the (1, 1) -radical inverse function defined by (fT2l) on N 1,1 is 

M 

13 



r 



^\X 




-^^^ 


" 



r 



Figure 2: The functions i//o and i//i associated to {<^f j} 

with the same coefficients akiji) of the binary representation n given by (fT3l) for 
Z7 = 2. 

By Theorem [53l the Kakutani-Fibonacci sequence of points {i^f j} coincides 
with the sequence {0i i(n)}. 

The following table shows the construction of the first 12 points of {<^"i}, 
where the first column contains the first 12 elements of Ni i. 
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0^ 


[0]2 = 


= -^ 


0.0 


= ■<^l,l(0)]2 


-> 





= ^/,i 


1^ 


:i]2 = 


= 1 -^ 


0.1 


= <^l,l(l)]2 


— )■ 


r 


= ^?i 


2^ 


2] 2 = 


= 10 ^ 


0.01 


= [<^m(2)]2 


-)■ 


f 


= ^?i 


4^ 


4] 2 = 


= 100 ^ 


0.001 = 


= [<^i,i(4)]2 


— )■ 


f 


= ^{a 


5^ 


5]2 = 


= 101 -^ 


0.101 = 


= >i,i(5)]2 


— )■ 


Y+Y^ 


= ^L 


8^ 


■8]2 = 


-- 1000 ^ 


0.0001 = 


= >i,i(8)]2 


-> 


/ 


= ^f.i 


9^ 


9]2 = 


= 1001 ^ 


0.1001 = 


= [<^m(9)]2 


— )■ 


7+/ 


= ^Zi 


10^ 


:io]2 = 


= 1010 ^ 


0.0101 = 


= [<^i,i(10)]2 


— )■ 


f + / 


~ ''1.1 


16^ 


;i6]2 = 


= 10000 -> 


0.00001 = 


= [0i,i(16)]2 


-)■ 


f 


~ ''1.1 


17^ 


:i7]2 = 


= 10001 -^ 


0.10001 = 


= [<^i,i(lV)]2 


— )■ 


Y+f 


_ ;:ro 
- '=1.1 


18-> 


[18]2 = 


-- 10010-^ 


0.01001 = 


= [<^i,i(18)]2 


— )■ 


Y^ + f 


_ ell 
~ '=1.1 


20^ 


20] 2 = 


= 10100 -^ 


0.00101 = 


= [<^ia(20)]2 


— )■ 


Y^j^f 


_£1'2 
~ '=1.1 
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